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1 .  INTRODUCTION 


This  report  contains  a  simplified  calculation  of  the  absorption  and  photo¬ 
luminescence  spectrum  of  modulation-doped  gallium  arsenide  for  energies  below 
this  material's  direct  band  gap.  The  calculation  was  motivated  by  the  results 
of  optical  experiments  [1]*  on  these  systems  which,  in  some  cases,  revealed  a 
two  peak  structure  in  the  photoluminescence  spectrum  whose  origin  was  unclear. 
The  work  described  below,  representing  a  first  attempt  at  clarifying  this 
situation,  is  divided  into  three  sections.  Section  2  is  an  overview  of  the 
problem  beginning  with  background  material  and  ending  with  a  description  of 
the  electronic  states  used  in  the  calculation.  In  section  3.  the  absorption 
spectrum,  the  photoluminescence  spectrum,  and  the  relationship  between  them 
are  derived.  Finally,  section  4  contains  a  discussion  of  the  theoretical  pre¬ 
dictions  along  with  critical  comments  on  various  simplifications  made  in  the 
analysis. 

2.  OVERVIEW 

The  systems  from  which  photoluminescence  spectra  were  obtained  consist  of 
approximately  l-pm-thick  wafers  of  gallium  arsenide  which  were  grown  by  mole¬ 
cular  beam  epitaxy.  During  growth  these  samples  were  selectively  doped  with 
donor  (Si)  and  acceptor  (Be)  impurities  at  densities  in  the  neighborhood  of 
10*®  cm~®.  A  typical  doping  profile  appears  in  figure  1,  where  a  positive 
ordinate  refers  to  donor  (n)  doping  and  a  negative  one  refers  to  acceptor  (p) 
doping.  The  abscissa  denotes  position  in  the  sample  measured  along  a  cubic 
axis  (GaAs  is  a  cubic  crystal  with  the  zinc  blende  structure).  The  periodic 
nature  of  the  doping  profile  has  led  some  to  refer  to  these  systems  as  super- 
lattices  because  of  the  additional  periodic  potential imposed  on  electrons  by 
the  ionized  donor  and  acceptor  atoms. 


DOPING 


Figure  1.  Doping  profile  along  (1,0,0)  direction 
of  gallium  arsenide  crystal. 


,  ‘References  appear  at  end  of  text. 

^Strictly  speaking,  the  potential  Is  not  periodic.  Since  the  Impurity  atoms  are  thought  to 
occupy  atomic  sites  in  a  random  manner  (at  least  at  these  doping  ^.ensltles),  the  added  potential 
can  only  be  viewed  as  periodic  in  the  growth  direction  If  one  Is  willing  to  average  over  the 
random  distribution  of  Impurities  In  the  transverse  directions.  This  point  will  be  addressed 
later. 

I 
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In  order  to  understand  the  optical  properties  of  this  system,  one  must 
have  a  clear  picture  of  the  electronic  states,  and  for  this,  a  brief  qualita¬ 
tive  description  of  the  effects  of  doping  in  a  single  crystal  of  GaAs  is 
appropriate.  Consider  first  the  case  of  a  single  donor.  Such  an  impurity 
differs  from  the  atom  it  substitutes.  In  particular,  when  neutral,  it  has  one 
extra  electron  occupying  a  localized  orbital.  This  electron  can  be  weakly 
bound  to  the  impurity:  its  energy  might  be  only  6  meV  below  the  bottom  of  the 
conduction  band.  In  such  a  circumstance,  the  effective  mass  theory  [2]  tells 
us  that  all  characteristics  of  the  localized  state  are  determined  by  the 
nature  of  the  low-energy  conduction-band  states.  For  instance,  the  "Bohr 
radius"  of  the  impurity  wave  function  is  given  by  ag  =  =  100  A,  where 

K  (=12.5)  is  the  static  dielectric  constant  and  m^  is  the  conduction-band 
effective  mass  (m^  =  0.067  m;  m  =  electronic  mass)  of  gallium  arsenide.  At  a 
doping  level  of  10‘*  cm“\  the  impurities  are  separated  by  an  average  distance 
of  about  10“®  cm  (-100  A),  so  that  the  wavef unctions  of  impurity  states  on 
nearby  sites  have  a  substantial  overlap.  This  means  that  electrons  are  no 
longer  localized  to  one  impurity  but  are  able  to  move  about  throughout  the 
crystal,  much  like  normal  conduction-band  electrons.  This  mobility  enables 
electrons  to  more  effectively  screen  the  impurity/ electron  interaction,  weak¬ 
ening  its  effect.  The  density  of  states  (DOS)  in  such  a  heavily  doped  system, 
with  the  DOS  in  the  undoped  case,  is  given  schematically  in  figure 
The  heavy  doping  changes  the  behavior  of  p(e)  near  the  band  bottom 
in  the  undoped  case  p(e)  behaves  as  [e  -  6^(0)  while  in  the  doped 

p(e)  gets  a  Lorentzian-like  contribution  to  the  intrinsic  behavior, 
decreasing  the  effective  band  minimum  by  a  small  amount  (a  few  milli-eleotron 
volts)  and  adding  to  the  overall  DOS  above  the  minimum.  As  is  evident  from 
figure  2(b),  the  situation  is  altogether  different  at  these  doping  densities 
for  acceptor  impurities.  Unlike  donors,  acceptors  are  neutral  when  their 
localized  state  (which  lies  just  above  the  valence  band  edge)  is  unoccupied. 
Theory  predicts  the  acceptor  state  in  GaAs  to  have  a  Bohr  radius  of  about  18  A 
(my  ;  0.35  mQ),  leading  to  very  little  overlap  between  acceptors  separated  by 
100  A.  The  states  would  retain  their  localized  nature  and  would  contribute  a 
relatively  sharp  peak  to  the  DOS  function  at  an  energy  above  the  valence  band 
edge  equal  to  their  "binding  energy." 


along 

2(a). 

e,(0): 

case 


Having  discussed  these  simpler  cases,  we  are  now  in  a  position  to  consider 
a  p-n  junction  and  then  a  superlattice.  Figure  3  shows  the  energy  bands 
versus  position  along  a  (1,0,0)  direction  in  GaAs  which  is  p-doped  to  the  left 
and  n-doped  to  the  right  of  z  =  0.  The  horizontal  dash-dot  line  represents 
the  Fermi  energy  of  the  system  whereas  the  heavy  dashed  line  represents  the 
acceptor  impurity  levels.  The  manner  in  which  the  charges  redistribute  them¬ 
selves  is  responsible  for  the  appearance  of  figure  3-  When  the  p-  and  n-doped 
systems  are  joined  together  at  z  =  0,  the  mobile  electrons  in  the  extended 
impurity  states  (z  >  0)  move  into  the  p-type  region  (z  <  0)  occupying  acceptor 
levels.  By  doing  this,  electrons  can  initially  lower  their  energy  by  -1.5  eV 
(the  energy  gap  Eg  in  GaAs).  As  more  electrons  fill  states,  the  vicinity  of 
the  junction  changes  its  character:  the  n-type  region  becomes  positively 
charged  and  the  p-type  region  becomes  negatively  charged.  This  charge  separa¬ 
tion  creates  an  electric  field  which  eventually  prohibits  further  charge 
flow.  The  total  depletion  region  (where  the  bands  "bend")  in  this  system  is 
about  650  A  wide.  This  width  leads  to  electric  field  strengths  of  about 
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Figure  2.  (a)  Conduction- 

band  density  of  states  in 
doped  and  heavily  doped 
n-type  gallium  arsenide 
and  (b)  valence-band  den¬ 
sity  of  states  in  doped 
and  undoped  p-type  gallium 
arsenide . 
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Figure  3.  Energy  band  diagram 
of  gallium  arsenide  p-n  Junc¬ 
tion  at  1  K  10“  cm“*  doping 
density. 
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1.5  V/650  A  -  10®  V/cm.  While  these  fields  are  large,  they  are  not  large 
enough  to  invalidate  the  semi  classical  theory  [3]  of  electron  dynamics  (which 
goes  into  the  arguments  leading  to  fig.  3)-  This  theory  allows  us  to  assume 
that  the  spectrum  of  single-particle  states  at  a  position  is  the  same  as 
that  found  at  Z2,  with  the  exception  of  an  overall  constant  shift  in  all 
energies  by  eC(t)(z2)  -  4i(zi)]  where  e  is  the  magnitude  of  the  electronic  charge 
and  <()  is  the  electrostatic  potential  set  up  by  the  ionized  impurities.  The 
validity  of  the  semiclassical  approach  implies,  among  other  things,  that  the 
DOS  above  the  conduct  ion- band  minimum  in  figure  3  is  like  that  shown  (roughly) 
in  figure  2(a)  for  all  values  of  z.  In  particular,  a  continuum  of  single¬ 
particle  states  begins  at  the  band  bottom  even  for  a  z-coordinate  within  the 
depletion  region. 

We  are  now  in  a  position  to  consider  the  nature  of  the  electronic  states 
in  a  crystal  of  GaAs  which  has  been  doped  according  to  the  profile  in  figure 
1.  Because  the  n  and  p  doping  densities  are  equal  and  because  the  doping 
period  5,  is  significantly  less  than  the  normal  depletion  width  (650  A)  in  a 
GaAs  p-n  junction  at  these  doping  densities,  the  lowest  energy  configuration 
of  this  system  would  correspond  to  completely  ionized  donors  and  acceptors: 
no  region  of  the  superlattice  will  be  charge  free.  The  results  of  a  semi¬ 
classical  approach  to  obtaining  the  electronic  structure  for  this  system  are 
displayed  in  figure  A.  The  conduction-band  bottom  (and  valence-band  top)  in 
this  figure  is  obtained  by  adding  to  the  unperturbed  band  structure  e^0)(k) 
the  potential  energy  U(z)  =  This  potential  energy  can  be  obtained  by 

"smoothing  out"  the  localized  ions  to  form  alternating  layers  (of  thickness 
1/2)  of  uniform  positive  and  negative  charge  density.*  We  then  need  to  solve 
V^4i(z)  =  -4itp(z)  with  p(z)  =  ±eN^,  where  the  upper  (lower)  sign  refers  to  n- 
doped  (p-doped)  regions  and  is  the  donor  (and  acceptor)  density.  This 
calculation  is  simple,  and  leads  to  figure  4,  with  Uq  given  by  Uq  =  irN^e*)l^/8< 
(which  for  the  parameters  of  fig.  1  equals  290  meV). 

Unfortunately,  the  semiclassical  approach  is  easily  shown  to  be  inadequate 
in  this  case.  A  calculation  of  the  kinetic  energy  (KE)  required  by  an  elec¬ 
tron  localized  in  one  of  the  conduction-band  wells  of  figure  4  gives  KE  -  3 
meV.  This  energy  would  represent  the  typical  separation  between  low-lying 
states  inside  the  well,  leading  to  a  DOS  far  below  that  predicted  by  the  semi¬ 
classical  picture  where  the  states  are  more  dense  by  a  factor  of  at  least 
500.  We  are  thus  forced  to  employ  the  effective  mass  theory  in  order  to 
obtain  the  electronic  states.  This  has  been  done  for  the  conduction-band 
states  by  Ruden  and  Dohler  [4].  Their  results  show  that  near  the  bottom  of 
the  wells  in  the  conduction  band  (and  the  top  of  the  valence-band  wells),  the 
electronic  states  are  harmonic  oscillator-like  in  the  z-direction  and  free 
electron-like  in  the  transverse  directions.  An  electron  in  one  of  these 
states  is  confined  in  the  z-direction  and  free  in  the  x-y  plane.  The  effec¬ 
tive  oscillator  frequency  is  easily  obtained  via  the  identification 
1/2[mQa)QZ^)  =  1 /2[4iTN^je*/K]z^ ,  giving  coq  =  (ATTN^je^/mpK) '  .  The  valence  band 

*The  approximation  of  replacing  the  Impurity  potential  by  Its  average  probably 

a  good  one  because  the  typical  difference  between  the  two,  ~  <t>(r)|>,  Is  less  than  the 

kinetic  energy  gained  from  one  electron  In  Its  delocalization  (which  Is  proportional  to  the 
impurity-impurity  overlap). 


is  more  complicated  because  in  GaAs,  the  valence-band  maximum  (at  q  =  0)  is 
four-fold  degenerate  (Including  spin).  In  addition,  the  effective  masses  are 
quite  different  from  the  conduction-band  mass.  These  two  facts  complicate  the 
structure.  In  the  present  work  this  complication  is  ignored  (comments  regard¬ 
ing  its  importance  will  be  made  in  sect.  ^4).  In  figure  5  a  schematic  of  the 
electronic  structure  is  given  along  with  a  DOS  profile  in  the  conduction  and 
valence  bands  (the  valence-band  profile  is  included  only  to  emphasize  its 
qualitative  difference  from  the  DOS  profile  of  fig.  2(b)).  Note  in  particular 
the  step-like  nature  of  these  profiles,  which  is  characteristic  of  a  two- 
dimensional  system  (the  DOS  in  a  two-dimensional  system  of  free  electrons  is  a 
constant).  The  acceptor  state  energies  are  represented  by  a  dashed  line  whose 
position  should  not  be  taken  too  seriously  since  our  treatment  of  the  valence- 
band  states  [5]*  has  been  approximate. 


Figure  M.  Energy  band  diagram  of  system  in  semi- 
classical  approximation. 


Figure  5.  Energy  band  diagram 
of  system  in  effective  mass 
approximation  Including  den¬ 
sity  of  states  profiles  in  the 
conduction-  and  valence-band 
wells.  For  clarity,  finite 
width  of  mlnibands  has  been 
ignored. 


-  CONTINUUM - 


'///'’ 

VALENCE  BAND 


*A  preliminary  calculation  of  acceptor  binding  energies  versus  position  In  a  superlattice 
(which  Ignores  valence-band  structure)  was  made  by  Leavitt  and  Simpson  [5]  and  supports 
qualitatively  the  positional  dependence  of  the  binding  energies  which  Is  displayed  In  figure  5. 
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In  light  of  this  electronic;  structure,  the  novel  optical  'and  electrical 
properties  of  this  system  are  readily  understood  in  a  qualitative  way.  The 
ground  state  has  all  acceptor  levels  filled  with  electrons,  and  the 
conduction-band  states  are  empty.  If  one  excites  electron-hole  pairs  with  an 
above-direct-band-gap  cw  laser  (e.g.,  a  He-Ne  laser),  electrons  will  scatter 
off  phonons  on  a  10“*^  s  time  scale  [6],  eventually  relaxing  into  the  lowest 
conduction-band  states  possible.  Similarly,  holes  will  relax  to  the  highest 
acceptor  states  available.  The  radiative  recombination  time  for  electrons  and 
holes,  which  in  a  uniformly  doped  semiconductor  is  of  order  10“’  s  [6],  is 
dramatically  increased  because  of  the  reduced  overlap  between  low-energy 
electron  states  and  acceptor  states  whose  wave  functions  are  separated  by  up 
to  a  superlattice  half  period  (200  A  in  this  system).  As  a  result  of  this 
reduction  in  recombination  rates,  a  cw  laser  will  initially  create  electron- 
hole  pairs  more  rapidly  than  they  can  recombine,  leading  to  a  filling  of 
electron  states  in  the  conduction  band  and  hole  states  in  the  acceptor  levels. 
As  this  filling  proceeds,  the  recombination  rate  rapidly  increases  for  two 
reasons.  First,  higher  oscillator  states  in  the  well  extend  further  out, 
giving  them  greater  overlap  with  acceptor  wave  functions.  Second,  and  more 
important,  as  the  steady-state  nonequilibrium  charge  density  in  the  conduction 
band  (and  acceptor  levels)  increases,  the  effective  potential  seen  by  an 
electron  is  weakened  by  the  carrier's  increased  ability  to  screen  the  space 
charge  field  of  the  ionized  impurities.  Any  weakening  of  this  potential  will 
dramatically  increase  conduction-band-state/acceptor-state  overlap,  resulting 
in  higher  radiative  recombination  rates.  When  the  generation  and  recombina¬ 
tion  rates  are  equal,  a  steady  state  will  be  reached.  This  steady  state  will 
be  characterized  by  an  effective  potential  that  differs  markedly  from  the 
potential  seen  in  the  absence  of  excitation.  This  means  that  each  pump  laser 
intensity  will  correspond  to  a  specific  steady-state  charge  density  in  the 
conduction  band,  n^^^,  and  its  corresponding  effective  potential.  Figure  6 
schematically  shows  the  situation  one  obtains  in  this  system  at  two  different 
pump  laser  intensities  (Ij^  >  1^).  The  higher  charge  density  in  b  results  in  a 
substantially  weakened  potential  and  increased  effective  gap  Eg. 

(a) 


(b) 
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Figure  6.  Energy  band  dia¬ 
grams  representing  two 
steady-state  nonequilibrium 
configurations  of  super¬ 
lattice.  Pump  laser  Inten¬ 
sity  which  maintains  steady 
state  Is  higher  in  (b)  than 
In  (a),  resulting  in  a  larger 
effective  gap. 
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de.scribed  in  tne  effective  mass  theory  by  wave 
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Here  is  tne  risen  functior.  correspondi  ng  to  the  lowest  energy  conduc¬ 

tion-band  state  and  F.,  r,  ;3  a  modulating  function  that  satisfies  the  equation 


)-  JiL 
)  ^ 
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(  P) 


eF  (r) 
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(2) 


where  e  is  the  energy  measured  with  respect  to  the  conducti on-band  effective 
potential  minimum  (at  z  =  0}  and  m,..  is  the  conduction-band  effective  mass. 
Because  of  the  translational  invariance  of  this  Hamiltonian  in  directions 
normal  to  the  z-axis,  F^,  ■'an  be  written 


(3) 


where  k  =  ,  C with  k,  =  '2-n/L)n^  on^  =0,  ±1,  ±2,  i  =  x,y). 

Unless  otherwise'  stated,  k  will  always  refer  to  a  two-dimensional  vector 
(k^  =  0).  We  taK<>  a  system  composed  of  periods  of  length  H  such  that  ?,Nj_  = 
L.  The  cross-sectional  area  of  each  layer  is  ,  so  that  the  total  volume  of 
the  system  is  a  =  L’.  Substituting  equation  (3)  into  equation  (2)  leads  to 
the  harmor.i oscillator  equation  for  n 


(4) 


whose  solutions 


an  be  w ’’it  ten 


where  b  =  (•h/mgWQ)  ,  H^(y)  are  Hermite  polynomials,  and  v  labels  the  oscil¬ 

lator  energy  levels.  The  eigenvalues  e  in  equation  (^4)  take  the  form  e^(k)  = 
+  •ftaiQCv  +  1/2).  We  have  chosen  to  normalize  the  wave  function 
according  to  the  requirement 

J  drltj)(r)|  *  =  1 


^  /n  =  1  . 

where  N  is  the  number  of  'unit  cells  in  the  system  and  the  second  integral  is 
over  the  volume  occupied  by  a  unit  cell  Qq. 

A  particularly  simple  form  has  been  chosen  for  the  acceptor  states:  one 
which  ignores  entirely  the  previously  mentioned  complications  in  the  valence 
band  by  treating  the  highest  energy  valence-band  state  as  nondegenerate 
(except  for  spin).  The  wave  functions  are  of  the  form 

iligCr)  =  FY(r)Uy^o('^)  » 

where  u^  highest  energy  valence-band  Bloch  function,  and  Fy  is  a 

modulating  function*  satisfying  a  hydrogen-like  Schroedinger  equation.  For 
simplicity,  Fy  is  taken  to  have  the  hydrogenic  form 


where  a  is  the  acceptor-state  Bohr  radius  and  r  measures  position  with  respect 
to  the  impurity. 

3.  CALCULATION 

We  begin  an  analysis  of  this  system's  optical  properties  by  determining 
its  absorption  coefficient.  In  this  calculation  we  assume  that  a  monochroma¬ 
tic  plane  wave  is  incident  on  the  superlattice  which  is  in  its  ground  state. 
The  Hamiltonian  of  the  system  is  given  by  H  =  Hq  +  ,  where  Hq  gives  rise  to 

the  harmonic  oscillator-like  and  acceptor  states  described  earlier,  and  , 
given  by  =  (e/mc)A*P,  represents  the  interaction  of  an  electron  with  an 
external  electromagnetic  field.  The  field,  which  is  assumed  to  be  linearly 
polarized  in  the  e  direction  (g*z  =0)  and  propagating  through  the  sample  in 
the  z-direction,  can  be  expressed  in  the  form 


E  =  Eo(a))c  +  c.c.  , 


*More  properly,  one  would  have  ♦,(r)  =  I  Fy(r)u  «  (r)  for  the  acceptor  level,  where  the  sum 

u  ’  • 

on  (I  is  over  all  degenerate  (q  =  0)  valence-band  states  (and  maybe  even  other  nearby  bands)  and 
the  Fy(r)  are  modulating  functions  satisfying  a  set  of  coupled  differential  equations.  See 
reference  2. 


where  q  =  ((D/c)n  and  n  is  the  refractive  index  of  GaAs.  This  field  is  deriv¬ 
able  from  the  vector  potential 


A  =  Ao(a))ee^  +  c.c.  ,  (9) 


I  with  E  =  -(1/c)  3A/3t  and  AgCu)  =  oEQ(a))/iaj. 

I  For  absorption  we  are  interested  in  the  rate  at  which  induces  transi¬ 

tions  between  the  acceptor  levels  and  the  conduction-band  states.  Labelling 
an  acceptor  state,  which  is  centered  at  (0,  0,  Zq) ,  by  "a"  and  conduction-band 
state  by  "\;,k"  we  obtain,  using  Fermi’s  Golden  Rule  [7],  the  rate  of  induced 
transitions  from  |a>  to  [v.lO; 


(10) 


where  e^Czg)  is  the  energy  of  the  acceptor  state  (measured  with  respect  to  the 
effective  potential  minimum  of  the  conduction  band).  Since  we  want  to  deter¬ 
mine  the  rate  at  which  energy  is  taken  out  of  the  radiation  field,  we  should 

multiply  this  by  fiw  and  then  sum  over  all  final  states  |v,k>.  The  delta 

function  will  ensure  that  only  states  with  the  proper  energy  are  included  in 

the  sum.  (Inspection  of  fig.  7  shows  that  for  a  given  tiw,  a  set  of  pairs 

(v,k^}  will  conserve  energy.)  Therefore  we  can  write 

W  =  tioo  I  I  ^  l<v,k|H^(u,)|a>l*  6(e^  ^  -  e^(ZQ^  -  fiuj)  ,  (11) 


(  f 


Figure  7.  Conduction-band-well  effective  oscillator 
and  density  of  states  along  with  a  single  hydrogenlc 
acceptor  wavefunctlon  located  a  distance  Zq  from  well 
minimum. 


where  W  is  the  rate  at  which  energy  is  absorbed  from  the  radiation  field  via 
transitions  out  of  one  acceptor  level.  If  we  use  the  relation 
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de  ,  ‘'v,k 
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L*m 


and  include  a  factor  of  2  for  spin,  equation  (11)  becomes 


2L*m  0) 
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where  k* 

0  for  x'^< 


=  (2m^/-h=')  jhu)  -  fi(*)o(v  +  1/2) 

0. 


"  "a(=^0)( 


and  0(x) 


(12) 


for  X  >  0  and 


In  order  to  determine  the  effects  of  all  the  impurity  states,  we  can  refer 
to  figure  8,  where  we  have  schematically  plotted  the  acceptor  levels  versus 
position  along  with  the  neighboring  conduction-band  wells.  The  contribution 
from  transitions  out  of  an  acceptor  to  its  neighboring  two  wells  will  be  w(zq) 
+  w[8,  -  Zq].  The  number  of  acceptors  in  the  slab  with  thickness  dZQ  at  Zq  is 
L^dZQN^.  We  can  then  integrate  over  Zq  from  i/4  to  3i/4  and  multiply  the 
result  by  Nj^,  the  number  of  periods,  to  get 


,3A/4 


"t  -  =^o)] 


Figure  8,  Acceptor  energy  profile  along  with  neigh 
boring  conduction-band  wells. 


Because  of  the  obvious  symmetry  about  z  =  1/2,  this  becomes 
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■m  , 


1/4 


2N 


iV'  “(I  *  *)  • 


Using  =  L  and  dividing  by  L*,  we  obtain  the  total  rate  of  energy  loss  per 
unit  volume 


^otal  ^/4  a  , 

-TT-  -  —  ^2  "  V  • 


(13) 


This  rate  is  simply  connected  to  the  absorption  coefficient  aCu).  Imagine  a 
wave  entering  a  slab  of  thickness  t.  On  entering,  the  wave's  intensity  is  Iq 
(  =  (nc/8iT)  |Eq|*),  while  on  leaving  it  is  lQe~“'' 

-  -  ....  .  _  _  «  _  “Cttf  *  —  — 


So,  the  net  energy  into  the 


slab  per  unit  time  is  IqL^  -  Ige  “  L*  -  IgatL*.  This  gives  Iga  for  the  net 
energy  into  the  slab  per  unit  time,  per  unit  volume.  Equating  this  to  the 
expression  above  gives 


2N  i/^ 

=  TT  J  dx  W(|  +  x)  . 
0^  -Jl/4  '  ' 


(14) 


To  proceed  further,  we  need  to  evaluate  W (Zq)  as  given  in  equation  (12). 
Consider  first  the  matrix  element  <v,k^|H^ (to) | a>.  This  can  be  written 

<v,kjH^(a.)|a>  =  J  dr  u;^Q(r)F*(r)(^  A(u,) -p) u^^ g(r)Fy (r  -  z^z)  .  (15) 


If  we  now  make  the  replacements 

r  ^  ^  P  .  Jo  ^  I  Jp  ' 

1  W  . 

and,  in  the  spirit  of  the  effective  mass  theory,  ignore  changes  in  the  mod¬ 
ulating  functions  on  the  length  scale  of  a  lattice  constant,  we  obtain 


<v,k  |H  (o))  I  a>  =  — 
V  1  ‘  me 


A(„).jl  F*(R.)PF,(r^  -  /  dp  u;,o(P)\.o(P) 


(16) 
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The  first  term  in  braces  is  zero  because  in  GaAs,  the  functions  u^  q  and  u^  q 
have  different  symmetry*  and  so  are  orthogonal.  Keeping  only  the  second  term 
and  using  AqCw)  =  cEQ(a))/ia)  leads  to 

eE-  \  1  c 

/^dRF^*(R)F,(B-z„z)  .  (17) 

where  we  have  defined^  Po,v^^^  “  ^  /ilO  integral  in 

equation  (17)  can  be  done ’more  easily  if  we  first  Fourier  transform  Fy,  re¬ 
placing  it  with 

Fy(x)  =  I  C(q)e^‘***  . 

q 

The  transform  C(q)  is  easily  shown  to  be 

r,  ^  1 

^  \  n  /  (1  +  ’ 


and  when  used  in  equation  (17)  gives 
1  >  ®Eo  , 


<v,ky  |Hi  (oj)  |a> 


(e«Pc,v(0))  I  C(ky,qz)nv(-qz)e 
Qz 


■i  Qz^O 


where 


-iq  z 


1  f”  ^z 

1  1  -  r  dz  n  (z)e  and  C 
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Consider  now  the  calculation  of  hy(qz)’  define 


N  =  /l  2"v  tt”*/-  b"*/* 

V 


and  use  the  dimensionless  variable  y  =  z/b,  we  can  write 


(%)  ■  -r  L  “y 


*s-11ke  and  p-llke,  respectively. 

tBecause  of  the  normalization  convention  chosen  for  the  u's,  the  factor  of  N  Is  needed  to  make 
the  result  finite. 
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Although  this  integral  can  be  done,  it  is  more  useful  at  this  point  to  define 
hv(Qz)  terms  of  derivatives  of  the  Hermite  polynomial  generating  function. 
For  this  we  define 


I  (a)  =  f  dy 

\1  *  -09 


■y*/2  -iay 


e  H  (y) 

V 


and  its  generating  function  F(a,a): 

»  I  (ci)s'^ 

F(3,0.)  .  I  -^4;, - /  dy 
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Since  \  — — j —  =  e  ^  ♦2sy  obtain  after  doing  the  integral, 


F(a,a)  =  e“  .  (22) 

Using  this  expression  along  with  the  generating  function  allows  for  a  deter¬ 
mination  of  Iy(a)  with  the  rule 


I  (a)  =  -  F(s,a)  ^ 

V  s=0 

Combining  these  results  in  equation  (20)  then  gives 
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After  n  (“Q  )  is  substituted  into  equation  (19),  the  matrix  element 
becomes  ^ 


<v,k^|H^  (o))  la> 


=  (e-P  (0))  N  —  (e^'f(3))  .  (25) 

ioim  V  c,v  /  /^  V  V  /  s=0 


where 


f(s)  =  f  dv 

^  —OB 


e  ''  ^^co3(v(23  + 


This  resul^  can  be  further  simplified  by  using  the  chain  rule  to  replace 
d'^/ds'^  (e®  f(s))|g^Q  by 

y  _ Vj _  d _  f*  f  Q  I 

J  (v  -  2p)!p!  ^jgV-2p  3-0  • 
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We  then  obtain 


<v;.k  |H  (ai)|a>  =  (e-P  (0))  N  I  - y-;.— ^  f(s) 

v'  1  '  luim  \  c,v  /  y—  V  p  (v  -  2p)!p!  ^gV-2p 

Also,  if  we  let  w  =  (a/b)v  and  u*  =  1  +  fg  can  be  written 

OH  H  ,  -(b*/2a*)w*  . 

f(3)  =  -  ^  r  dw  ^ - C03(rw) 

a  dy*  ■' 
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(28) 


differentiated  with  respect  to  y*  to  obtain 


where  f  »  (2b/a)s  +  Zg/a.  This  integral  can  be  done  [8]  and  the  expression 

to  y *  to  obtai n 
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where  X^^  =  [(b/a)y^  +  ZQ/b]//2  and  Erfc  is  the  complementary  error  func¬ 
tion  [8]."  We  can  now  use  these  results  to  obtain  for  W 


64a ’m  e* I E^| * 
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1  —  ©(xM 

,=0  \  v/ 


y  ^ 

v-2p 

p  (v  -  2p)!p 

1  O  »>v  *  '  ®  ' 

'  ds^'^P 

3»0_ 

(30) 


This  is  our  final  result  for  W,  which  can  be  used  in  equation  (14)  to  obtain 
the  absorption  coefficient  a((ii). 

Before  examining  these  results  in  any  greater  detail,  we  now  consider  the 
case  of  photoluminescence.  For  this,  we  assume  that  a  steady-state  non¬ 
equilibrium  condition  is  maintained  in  the  superlattice  which  produces  a  sharp 
quasi-Fermi  level  in  the  conduction  band.  Below  this  level,  all  conduction 
band  states  are  occupied.  Spontaneous  emission  of  photons  will  result  from 
transitions  of  these  electrons  into  empty  acceptor  states.  For  photolumines¬ 
cence,  we  need  to  determine  the  rate  at  which  these  photons  are  emitted. 
Detailed  balance  arguments  given  in  appendix  A  lead  to  a  relationship  between 
microscopic  induced  and  spontaneous  transition  probabilities,  which  can  be 
expressed  through  the  equation* 


dN 


dt 


*In  this  equation,  statistical  factors  are  assumed  equal  to  values  which  would  maximize  each 
ratio  (see  app  A) . 
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where  Nq  is  the  number  of  photons  in  mode  q,  \  of  the  radiation  field.  To 
obtain  the  photoluminescence  spectrum  R(e),  which  was  defined  earlier,  we  need 
to  multiply  equation  (31)  by  fioj,  sum  over  all  photon  modes  q,i  consistent  with 
their  energy  being  in  the  interval  de  at  e,  and  divide  the  result  by  the 
volume  of  the  system.  This  gives 
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where  the  prime  on  the  sum  is  a  reminder  of  the  energy  restriction.  But  Nq^x 
-  (u  ^^Q/'TUo),  where  is  the  q,X  contribution  to  the  energy  density  in  the 

fiela.’  If  we  use  this ’together  with  equation  (31).  we  obtain 
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The  quantity  “fiwdNq^  ;^/dt  |  ^  ^d  .abs .  rate  at  which  energy  is  absorbed  from 

the  radiation  fi^d,  which  is  precisely  the  quantity  lQa(m)fl  calculated 
earlier  (provided  we  take  q  =  (0,0, q)  and  choose  the  same  polarization 
e).  Thus  we  can  write 


I-a  .(a))fl 
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Using  this  result  in  equation  (33)  along  with  the  relation  u  =  nlQ/c  gives 
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The  only  dependence  of  aq  on  q  and  A  comes  from  the  square  of  the  matrix 
element  which  contains  ^'e  factor  |e;^(q)*P^;  y(0)|*.  Here  '  ^*2)  are 
the  orthonormal  polarization  basis  vectors  which  are  orthogonal  to  q  and  de¬ 
fined  according  to  some  arbitrary  (but  definite)  prescription.  With  this  in 
mind  we  consider  the  sum 
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Using  e  =  'fico  =  f\cq/n,  this  becomes 
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where  we  have  defined  tt 
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1  j  dq  €j^(q)ej^(q) . 


(37) 


One  can  replace  dq  by  de  and  easily  obtain  =  (l /3)qiqj6i j .  If  this  is 
used  in  equation  (37).  one  obtains 


2  to^n^de  1 


(38) 


Substituting  this  result  into  equation  (35)  then  gives 


R(e)  =  i(a))>  , 

^  c  q,A 


(39) 


where  <aq,x(‘*>)^  [=1/^ir  j  dflq  Oq^x^***))  the  angular  average  of  x^***^'  which 
is  given  by  equation  (1^)  for  the  case  q  =  (0,0,q)  and  polarization  =  e.  As 
shown  above,  <a(a))>  can  be  obtained  simply  from  expressions  for  a(u)  by  re¬ 
placing  the  quantity  |e*Pc,v(0)l^  which  appears  in  the  squared  matrix  element 
by  the  expression  1/3 


Equations  (14),  (30),  and  (39)*  are  the  central  results  of  this  work. 

Under  appropriate  circumstances  they  can  be  used  to  gain  insight  into  the 
photoluminescence  and  absorption  spectra  of  doping  superlattices. 


4.  PREDICTIONS  AND  DISCUSSION 

Before  the  results  obtained  above  are  examined,  a  few  comments  will  be 
made  regarding  their  applicability. 

The  absorption  spectrum  given  by  equation  (14)  represents  contributions  to 
absorption  from  electronic  transitions  between  localized  acceptor  states  and 
extended  (in  two  dimensions)  conduction-band  states.  When  liu)  is  large  enough 
to  induce  transitions  between  extended  valence-band  and  extended  conduction- 
band  states,  equation  (14)  will  no  longer  apply,  as  these  contributions  were 
ignored  in  the  analysis.  It  is  also  important  to  remember  that  this  formula 
is  applicable  only  when  the  exciting  field  is  reasonably  weak."'"  An  intense 
external  field  will  generate  electron-hole  pairs  more  rapidly  than  they  can 
recombine,  leading  to  a  finite  steady-state  conduct  ion- band  charge  density 
n^2)  an(j  altered  effective  potential.  In  this  instance,  the  potential  used 
in  obtaining  the  approximate  eigenstates  will  be  inappropriate,  and  equation 
(14)  will  provide  only  a  qualitatively  correct  account  of  the  spectrum. 

*Care  should  be  exercised  In  using  equation  (39)  for  the  photoluminescence  spectrum.  In  our 
calculation  of  qi(w)  we  considered  contributions  from  all  acceptor  levels.  For  photoluminescence, 
all  acceptor  states  cannot  be  final  states  since  some  are  already  occupied.  Instead,  one  should 
limit  the  range  of  Integration  In  equation  (14)  to  the  distance  required  to  contain  a  number  of 
states  equal  to  the  total  number  of  occupied  states  In  the  conduction  band. 

^In  view  of  the  nonlinear  optical  properties  of  this  system,  this  constraint  could  be  severe. 


In  the  case  of  photoluminescenee,  similar  considerations  will  apply.  In 
particular,  high  pump  laser  intensities  will  produce  steady  states  whose  cor¬ 
responding  effective  potentials  are  substantially  different  from  the  effective 
potential  used  in  the  calculation  (which  corresponded  to  n^^)  =  0).  Despite 
these  limitations  of  the  theory,  the  more  serious  electronic  structure  calcu¬ 
lations  of  Ruden  and  Dohler  [4]  support  the  qualitative  picture  presented  in 
the  model,  at  least  for  values  of  n^^^  which  are  not  excessive  ( >50  percent  of 
the  available  acceptor  charge). 

In  what  follows  we  examine  the  photoluminescence  spectrum  predicted  by 
equation  (39)  for  a  superlattice  characterized  by  a  doping  density  of  1  x  10'® 
cm“®  and  a  period  of  800  A.  We  will  consider  a  range  of  charge  densities  n^^^ 
since  the  spectra  can  be  rather  sensitive  to  this  density.  In  order  to  apply 
equation  (39),  we  need  to  determine  two  things.  First,  we  need  a  reasonable 
representation  of  the  manner  in  which  the  acceptor  binding  energies  vary  with 
distance  from  the  center  of  the  p  regions.  Variational  calculations  [5] 
support  qualitatively  the  use  of  a  binding  energy  profile  similar  to  the  one 
displayed  schematically  in  figure  8.  The  binding  energies  are  largest  in  the 
center  of  the  p  regions  and  decrease  to  small  values  near  the  edge  of  the  p 
regions  where  the  space  charge  field  is  largest  (they  can  never  go  lower  than 
the  V  =  0  oscillator  level  in  the  valence  band  since  this  marks  the  beginning 
of  a  continuum  of  states).  To  represent  this  state  of  affairs,  the  binding 
energies  were  modelled  by  the  expression  ejj(z)  =  e"^z^,  where  z  measures 

distance  from  the  center  of  the  p  regions  and  the  parameters  Y  and  e^xnax 
chosen  judiciously.  Finally,  a  choice  must  be  made  regarding  the  way  in  which 
empty  acceptor  states  are  distributed  (among  the  occupied  acceptor  states) 
when  charge  is  elevated  into  the  conduction  band  in  order  to  represent  a 
specific  steady-state  configuration.  We  will  assume  that  the  empty  acceptor 
states  are  those  acceptors  which  occupy  a  region  of  finite  thickness  centered 
on  the  p  regions  (i.e.,  the  highest  energy  states).  The  thickness  of  these 
regions  is  determined  from  the  nonequilibrium  charge  density  in  the  conduction 
band.  For  instance,  if  50  percent  of  the  available  acceptor  charge  (electrons 
occupying  acceptor  states  in  the  ground  state)  were  elevated  into  the  conduc¬ 
tion  band,  then  one  half  the  thickness  of  the  p  layers  would  be  neutral  (200  A 
for  a  superlattice  period  of  800  A).  We  remind  the  reader  that  in  such  a  case 
the  binding  energy  profile  would  change  since  the  "background"  potential  would 
be  flattened  in  the  central  halves  of  the  p  regions.*  The  Gaussian  profile 
function  given  above  is  flexible  enough  to  account  for  this  effect  qualita¬ 
tively  since  decreasing  Y  (or  £5(^3^^  would  flatten  out  the  profile. 

For  the  parameters  we  will  consider,  the  quasi-Fermi  energy  Ep  will  always 
lie  below  the  v  =  2  oscillator  level.  In  this  case  equation  (30)  simplifies 
considerably  and  can  be  combined  with  equations  (14)  and  (39)  to  give,  after 
some  manipulation,  the  expression 

R(e)  =  D^*F(^)  ,  (40) 

‘Another  aspect  of  the  binding  energy  profile  which  we  have  Ignored  In  the  present  calculation 
Is  the  Inevitable  dependence  of  the  average  acceptor  wavefunctlon  radius  on  Its  binding  energy 
and  hence  position.  Our  present  lack  of  information  regarding  the  acceptor  wavefuncti ons 
necessitates  this. 
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where  D  is  a  constant,*  C  is  the  photon  energy  in  units  of  fiojQ,  and  F(5)  is 
given  by  the  integral 


F(o  -  p  dw  {iSs(K5)[i  -  e(kjf)]  .  1  i;G(k;)[,  -  e(k;^)]}  . 

«/  w 


(41  ) 


where  0  is  the  step  function  defined  below  equation  (12), 

=  [2/b^])c  -  (n  +  1/2)  +  ea(a))/-Ra)oj  , 

and 


is  k^  with  5  =  Sf  =  ep/'I^O 


The  functions  [l  -  0(k^]]  prevent  photoluminescence  contributions  from  states 
above  the  quasi-Fermi  level  while  the  functions  ©(k^)  prevent  contributions 
from  below  the  nth  oscillator  levels.  The  functions  Iq  and  are  given  by 


0  2MQa 
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2a^p^ 


< 


e  e  +  e  e 


^  b"p?/2a"  ^ 

•rrb  1  Tp 

-  e  e 

ap 


.b  -^^1 

+  -  e  e 

api 


'Erfo(x,J  (i  ,,  . 

EErc(x,  j[j|;  -  «)j 


(43) 


b%V2a^ 

-  (Tr/2)‘/^  ^  e  (1  -  b%?/a^) 

ap'  1 


rpi  -xp  -Tp^  -xp 

e  e  -  e  e 


-J  . 


where  X  =  [(b/a)p  ±  w]/y2  and  r  =  (b/a)w.  To  define  Ep  we  decide  what 


vl  ^  '  "V  ■■  ■  ■  'F" 

fraction  of  the  available  acceptor  charge  is  to  be  placed  in  the  conduction 

band.  The  maximum  available  density  n^^^  (per  unit  area,  per  layer)  is  simply 


max 


given  by  1/2,  while  Ep  is  given  by  the  solution  to  the  equation 

(mpirfi"^  (ep  -  fia)o/2)  =  n^^^ 


*D  1s  given  by  0  =  (1024a/3  n’^’n)  (N^a*  )  (e*/linc’ XNj/lmTi)  |  Pp  ^  y  ( 0)  |  *  where  o  (equal  to  e'^/Tic) 
Is  the  fine  structure  constant. 


22 


.>  .'p  .I- 


i 


r4 


C. 

9 


V 


n: 


i 


when  Ep  lies  between  the  v  =  0  and  v  =  1  oscillator  levels  or 


+  (2m^/Tm*)  (cp  -  S/^nojo)  =  n^^^ 


when  Ep  lies  between  the  v  =  1  and  v  =  2  oscillator  levels.  In  these 
expressions  m^/irTi*  is  the  DOS  (per  unit  area)  in  a  two-dimensional  system 
[3].  Finally,  the  limits  of  integration  in  equation  (42)  are  given  by 


/o  -  (2),  (2)\  I 

w  =  (2  +  n  /n  j  "iTir  ■ 
1,2  \  max/  4b 


In  figure  9  a  sequence  of  photoluminescence  spectra  is  displayed  for 
successively  higher  values  of  Each  curve  is  labelled  by  the  ratio 

n(2)/n(2)  Because  of  the  theory's  inability  to  accurately  predict  the 


n( 2) /p( 2) ^ 

max 


Because 


theory' s 


?ctra  is  displayed  for 
labelled  by  the  ratio 
accurately  predict  the 
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Figure  9.  Theoretically  pre¬ 
dicted  photoluminescence 
spectra  for  superlattice  with 
parameters  given  In  figure  1 
(also  ”  20  meV  and  1  - 

1 /T*  -  75  A;  see  text) . 
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effective  gap,  all  effective  gaps  were  set  arbitrarily  to  1 , 3  eV.  Hence  only 
differences  in  the  horizontal  axis  are  significant.  The  vertical  scales  are 
in  arbitrary  units.  The  curves  were  obtained  by  numerical  evaluation*^  of  F(C) 
with  the  use  of  the  parameters  =  20  meV  and  \  =  75  A.  The 

nature  of  the  curves  is  readily  understood  if  attention  is  paid  to  two  general 
considerations: 

(a)  The  square  of  the  matrix  element  of  decreases  exponentially  with 
separation  between  the  "centers"  of  the  two  wavef unctions  involved.  This 
implies  that  for  a  given  acceptor  position,  higher  oscillator  levels  will  have 
significantly  larger  squared  matrix  elements  since  their  wavef unctions  extend 
further  out  from  the  oscillator  center. 

(b)  For  a  given  acceptor  position  and  oscillator  level,  the  squared  matrix 
element  of  will  decrease  with  increasing  -ho)  since  this  larger  energy  will 
go  into  conduct i on- band  wavefunction  phase  oscillations,  resulting  in  smaller 
overlap  with  the  hydrogenic  orbital. 

Let  us  consider  in  detail  the  features  of  figure  9,  curve  a.  The  thresh¬ 
old  for  photoluminescence  begins  at  1 . 3  eV  and  rises  exponentially  to  a  maxi¬ 
mum  at  1.312  eV.  This  rise  occurs  because  at  threshold  only  states  near  the 
top  of  the  acceptor  profile  (and  hence  further  away  from  the  oscillator  well) 
are  permitted  to  contribute.  After  hw  has  increased  to  above  1.312  eV,  the 
entire  range  of  acceptor  energies  is  included,  and  a  further  increase  in  hco 
leads  to  a  decay  in  the  contri  but  i  ons  due  to  consideration  (b).  This  decay 
will  continue  until  the  threshold  condition  for  the  next  (v  =  1)  oscillator 
level  is  reached.  One  then  gets  an  exponentially  increasing  contribution 
w.ni  ch  will  persist  until  either  the  quasi-Fermi  energy  is  reached  or  Tiw  has 
increased  by  0.012  eV  (  l  .  31 2  eV  -  1.300  eV),  where  another  decay  will  begin. 
This  is  what  happens  as  one  goes  from  curve  a  to  d.  As  the  quasi-Fermi  level 
increises,  the  second  peak  grows  relative  to  the  first  until  (in  d)  the  second 
peak  completely  .'dwarfs  the  contributions  from  the  v  =  0  level,  leaving  only  an 
asymmetric  peaw.  The  behavior  in  figure  9  is  generic  and  can  be  reproduced 
for  a  variety  of  parameter  sets. 

Figure  10  displays  experimental  photoluminescence  results  of  Simpson  et  al 
'1]  on  a  superlattice  sample  characterized  by  parameters  similar  to  those 
consi dered  theoretically.  The  photoluminescence  intensity  is  in  arbitrary 
units.  when  an  attempt  is  made  to  compare  the  theoretical  predictions  with 
this  spectrum,  a  serious  difficulty  arises.  In  the  experiments,  a  pump  laser 
is  employed  whose  penetration  depth  into  the  superlattice  is  approximately 
1000  to  2000  A.  We  are  thus  faced  with  the  difficulty  of  choosing  values  of 
n  21  for  each  successive  conduction-band  well.  Each  well's  effective  gap  (and 
relevant  matrix  elements)  will  then  be  changed  by  different  amounts,  and  the 
total  1  umi  nesceni^e  will  be  the  sum  of  contributions  from  all  layers.  Such  a 
calculation  involves  too  many  unknown  parameters  and  would  be  unconvincing. 


PHOIDLUMINESCENCE  SPECTRUM 
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Figure  10.  Experimental  photolurai nescence  data  of 
Simpson  et  al  [1 ], 

Because  of  this  difficulty,  we  will  postpone  serious  comparisons  between 
theory  and  experiment  until  a  less  restrictive  theory  (or  more  compatible  ex¬ 
periment)  exists. 

In  concluding  this  report  we  list  below  a  number  of  criticisms  and  com¬ 
ments  relating  to  this  model. 

1.  We  have  considered  only  the  case  where  =  N^,  but  this  precise 

equality  is  experimentally  unlikely.  If  were  greater  than  the  poten¬ 
tial  in  the  p  regions  would  be  flattened  out  in  the  ground  state.  On  the 
other  hand,  if  were  greater  than  N^,  the  n  region  potential  would  flatten 
and  the  ground  state  would  be  characterized  by  a  nonzero  conduction-band 
charge  density  The  qualitative  features  of  the  theory  would  remain 

unchanged  in  these  cases. 

2.  We  have  ignored  details  of  the  valence  band  in  the  analysis.  It  is 

hard  to  see  how  qualitative  differences  would  arise  because  of  this  neglect. 

Although  the  effective  mass  theory  would  present  more  formidable  mathematical 
problems  in  obtaining  the  acceptor  state  wavefunctions ,  the  wavef unctions 
would  still  be  qualitatively  hydrogenic,  with  some  distortion  in  shape-- 

possibly  flattened  spheres  with  the  thin  dimension  normal  to  the  layers  (these 
wavefunctions  could  also  be  obtained  variationally  if  desired).  The  only 
aspect  of  this  result  which  could  be  qualitatively  important  is  that  if  the 
long  dimension  of  the  wavefunction  became  sizable,  overlap  with  similar  wave- 
functions  might  produce  extended  (in  two  dimensions)  acceptor  states.  This 
possibility  remains  to  be  adequately  addressed. 

3.  We  have  ignored  the  effects  of  randomness  in  the  analysis.  It  is 

clear  that  since  the  acceptors  are  at  random  positions  throughout  the 
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p-regions  (with  an  average  separation  of  100  A),  an  electron  in  the  immediate 
neighborhood  of  a  specific  acceptor  will  feel  a  "background  potential"  (the 
potential  due  to  all  other  changed  acceptors  and  donors)  which  differs  from 
one  acceptor  to  another  at  a  fixed  z-coordinate.  This  variance  in  background 
potential  was  ignored  in  the  analysis,  and  only  the  average  background  poten¬ 
tial  was  included  (which  varied  quadratically  with  distance  from  the  center  of 
the  p-region).  It  is  difficult  to  estimate  the  size  of  this  variance  since 
charge  carrier  screening  of  the  potential  is  probably  substantial.  (A  simple 
model  including  screening  is  currently  being  studied  to  give  some  idea  of  this 
variance  in  background  potential.  The  results  of  this  investigation  will  be 
presented  in  a  future  publication.)  In  any  event,  this  phenomenon  will  lead 
to  a  distribution  of  binding  energies  at  each  coordinate  Zq.  One  would  then 
have  to  integrate  over  this  impurity  band  for  each  Zq,  leading  to  obvious 
changes  in  the  photoluminescence  spectrum. 

4.  In  the  analysis  we  ignored  the  fact  that  decreasing  binding  energy 
will  mean  increasing  acceptor  state  extension.  A  more  accurate  treatment 
would  have  allowed  the  acceptor  Bohr  radius  a  to  be  a  function  of  z. 

5.  Homogeneous  line  broadening  was  ignored  in  the  analysis  but  would 
result  in  a  smoothing  out  of  sharp  spectral  features  in  figure  9. 

6.  The  reader  is  reminded  that  the  spectra  displayed  in  figure  9  are  very 
sensitive  to  the  distribution  of  occupied  acceptor  states.  In  particular,  the 
10*°  factor  in  figure  9(d)  comes  primarily  from  the  fact  that  higher  densities 
provide  empty  acceptor  states  closer  to  the  conduction-band  wells.  A  differ¬ 
ent  acceptor  state  distribution  would  result  in  an  entirely  different 
spectrum . 
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APPENDIX  A.— RELATIONSHIP  BETWEEN  INDUCED  AND  SPONTANEOUS 
TRANSITION  PROBABILITY 


APPENDIX  A 


The  principle  of  detailed  balance  permits  one  to  establish  a  relationship 
between  the  microscopic  transition  probabilities  for  induced  emission  and 
spontaneous  emission  in  a  physical  system  without  resorting  to  the  otherwise 
unnecessary  procedure  of  quantizing  the  electromagnetic  field. 

Consider  a  superlattice  in  thermal  equilibrium  with  the  radiation  field, 
and  focus  attention  on  two  specific  states  in  the  superlattice,  say  "1"  and 
"2"  with  energies  e-j  and  €2  (<ei).  respectively.  We  assume  that  electrons 
making  transitions  between  these  states  will  be  accompanied  by  emission  and 
absorption  of  photons  of  type  q,X  changing  the  state  of  the  radiation  field. 
If  we  let  Nq  ^  represent  the  number  of  photons  in  this  mode  of  the  field,  we 
can  write  roV  the  total  rate  at  which  N  changes  because  of  occupancy 
changes  in  states  1  and  2: 


dt 


dt 


dN  , 

dN  , 

q.X 

q.X 

•  dt 

..  dt 

ind.em. 

spont.em. 

ind .abs . 


The  induced  absorption  and  emission  rates  are  proportional  to  the  number  of 
photons  present  and  to  the  probability  that  the  appropriate  initial  states  are 
occupied  and  final  states  unoccupied: 


dN 


SUA 


dt 

dN  , 
dt 


ind.em.  ®%,X^20  ^l) 

ind. abs.  "  "®%,X^1  0  ”  ^2) 


The  fact  that  both  have  the  same  coefficient  B  follows  from  microscopic  re¬ 
versibility.  On  the  other  hand,  the  spontaneous  emission  rate  simply  depenas 
on  the  probability  that  state  2  is  occupied  and  state  1  unoccupied: 


dN 


q,  X 


dt 


spent .em. 


■  ''■’20 
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When  thermal  equilibrium  prevails,  dN  /dt  =  0  or 

q ,  X 


0  - 

1  +  BN 

-  P.) 

=  bn  ,p,  (: 

1 

) 

\  V 

'  q.X  2V 

1/ 

q.X  U 

which  gives 


raO  -  b)] 


But  in  equilibrium  U  ^  w)/kTj  ^  little  algebra  reduces  the 

bracketed  term  to  -  1,  which  is  precisely  the  inverse  of  the  Planck 
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distribution  Nq^x  =  (g^w/kT  _  i)-i^  Thus  A  =  B,  and  for  a  system  in  which  we 
know  from  first  principles  the  rate  ind.abs  •  write 


dN 


q,X 


dN 


q.X 


dt 


dt 


=  B  = 


ind.abs , 


spont .em. 


q,X 


where  the  statistical  factors  are  assumed  equal  to  P-j  =1,  P2  =  0  for  induced 
absorption  and  P^  =0,  P2  =  1  spontaneous  emission. 
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APPENDIX  B.— FORTRAN  CODES  FOR  PHOTOLUMINESCENCE  CALCULATION 
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Appendix  B 
integral  given 


lists  the  Fortran  code  used  in  the  n'jmerical  evaluation  of  the 
by  equation  (^^i)  in  the  main  body  of  this  report. 


. . .  CALCULATION  . . . . 

Cl  SCREEN 

C|  WHITE  SYSIOUTPUT  "COMPILING" 

Cl  FORT  PLSPEC 

C|  WRITE  SYSIOUTPUT  "LINKING" 

Cl  LINK  PLSPEC.XYPLOT.HDLILIB/L.CMLIB/L.DISSPLA/L.GH/L 
C 

C  . .  BEGIN  MAIN  PROGRAM  . 

C 

IMPLICIT  REAL'S  (A-H.O-Z) 

REAL*4  SENERGYI 1000) .SPLSPEC( 1000) ,YB( 1000) 

CHARACTER* 1  NANS 

D I MENS I ON  ENERGY ( 1 000 ) . PLSPEC ( 1 000 ) . I WORK (100), WORK ( 400 ) 
COMMON /CONSTANTS /X . A , B . P I . SP 1 2 . STWO . SEP . EGPEXP . EBMAX . SP , ALPHA 
COMMON/EXTRA  A'PL I M . GAMMA 
EXTERNAL  F.DERFC.EA 

DEFINE  FUNCTIONS 

YY(X.W)=2.*(A/B)"2*<X-0.5+(EA<W)/SEP)) 

AMU(X.W)=DSqRT( 1 .+YY(X.W) ) 

EB ( 2 ) =EBMAX • DEXP ( - A LPHA • 2 • • 2 ) 

ENTER  DOPING  PARAMETERS 

TYPE  '.ENTER  DOPING  DENSITY  IN  UNITS  OF  1.0X10(18)  CM-3 
ACCEPT  ',DD 
DD=l .0D18»D0 

TYPE  '.ENTER  SUPERLATTICE  PERIOD  IN  UNITS  OF  ANGSTROMS' 
ACCEPT  '.SP 
SP=1 .OD-8'SP 

DEFINE  VARIOUS  CONSTANTS 

EG= 1 , 424 
AK=12.5 
AM0=9. 1005D-28 
AMV=0.50D0'AM0 
AMC=0.067D0'AM0 
HBAR=1 .05459D-27 
EC=4.80324D-10 
C0NV=1 .60219D-12 
A= ( ( HBAR/EC ) "  2 ) ' AK/AMV 
BBEEV=(  AMV/(  AMO' AK"2 ) ) '  1 3 . 60568D0 
BBMEV= 1 000 . DO'BBEEV 
BBE=CONV'BBEEV 
PI=4.D0'DATAN( 1  DO) 

SP I 2=DSqRT ( P I /2 . DO ) 

STW0=SQRT(2.D0) 

CALCULATE  OMEGAO,  OSCILLATOR  ENERGY  HBAR'OMEGAO, 

AND  ZERO  POINT  AMPLITUDE  B 

OMEGAO=EC'DSQRT( 4 . DO'PI 'DD )/DSqRT( AMC'AK ) 

SEP=HBAR' OMEGAO 
SEPEV=SEP/CONV 
SEPMEV= 1 000 . DO ' SEPEV 
TYPE  I .SEPMEV 
I  F0RMAT(/.20X, 'HBAR'OMEGAO=' .F7 . 3 . 'mev ' ./) 

B=SQRT( HBAR/( AMC'OMEGAO ) ) 

BANC=1 .008'B 
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C 

r  CALCULATE  UO  AND  EEG(EFF  ENERGY  GAP) 

C 

UO= ( AMC* ( OMEGAO*SP ) • *2 )/( 32 . DO*CONV ) 

EEG=EG-2.D0*U0 
IF(EEG.GE.O.DO)GO  TO  2 

TYPE  •.•STOP  IN  PLSPEC:  EFF  ENERGY  GAP  IS  NEGATIVE' 

STOP 

2  CONTINUE 

C 

C  DETERMINE  2-D  DENSITY  OF  STATES 

C  (NOTE  THAT  BECAUSE  OF  UNITS.  HEAR* *2  IS  TOO  SMALL 

C  AND  SO  HAS  BEEN  REPLACED  WITH  HEAR.  CORRECTIONS  ARE 

C  MADE  IN  APPROPRIATE  PLACES) 

C 

RH02D=AMC/( P I ‘HBAR ) 

ST ATES=RH02D • SEP/HBAR 
C 

C  DETERMINE  NON-EQUILIBRIUM  CHARGE  DENSITY* 

C 

AVAILC=DD*SP/2.D0 
TYPE  4.AVAILC 

4  FORMAT(/.'  MAX.  NON-EQUI LIBRIUM  CHARGE  DENSITY= ' .D8 . 3 .' cm-2 '. /) 

TYPE  •.•ENTER  ACTXIAL  NON-EQUILIBRIUM  CHARGE  DENSITY  AS  A  FRACTIO 
CN  OF  THE  ABOVE  STATED  MAXIMUM  (NUMBER  BETWEEN  0.  AND  I.)' 

ACCEPT  • . FRAC 
CHDENS=FRAC^AVA 1 LC 
C 

C  CONSTRUCT  MODEL  BINDING  ENERGY  PROFILE 

C 

TYPE  23.BBMEV 

23  FORMAT( •  ACCEPTOR  BINDING  ENERGY  IN  UNDOPED  CRYSTAL  IS* .F7.3. ’Be 
Cv*  ./.•  ENTER  MAX  BINDING  ENERGY  IN  SUPERLATTICE  IN  Bev'./) 

ACCEPT  •.EBMAX 
EBMAX=EBMAX •CONV/1 000 . DO 
TYPE  24. BANG 

24  FORMATC  OSCILLATOR  ZERO  POINT  AMPLITUDE  I S ' ,F7 . 2 . * angs t roms * . / . 
C*  ENTER  LENGTH  SCALE  FOR  BINDING  ENERGY  DECREASE  IN  angatroBS * ./) 

ACCEPT  •.ALPHA 
ALPHA=ALPHA/BANG 
ALPHA=1 .D0/( ALPHAS •Z) 

C 

C  DETERMINE  INTEGRATION  LIMITS  AND  MAX.  ACCEPTS  ENERGY 

C  RANGE  TO  BE  ENCOUNTERED 

C 

ACCTHK=FRAC^SP/2 . DO 

WM I N=0 . 5D0^ ( SP-ACCTHK )/B 

WMAX=0 . 5D0^ ( SP+ACCTHK )/B 

TYPE  • . •WMIN=' .WMIN. •  WMAX=- .WMAX 

WZ=WMAX-0 . 5DO^SP/B 

ERANGE=EBMAX-EB( WZ ) 

ERANGEME V= 1 000 . DO • ERANGE/CONV 
TYPE  333 . ERANGEMEV 

333  FORMAT ( ’  ENERGY  RANGE= ' , F6 . 2 , ' MEV ' ) 

C 

C  DETERMINE  FERMI  LEVEL 

C 

CHMAX=3 . DO •STATES 
I F ( CHDENS . GE . CHMAX ) THEN 

TYPE  •.STOP  IN  PLSPEC:  PROGRAM  WONT  PRESENTLY  HANDLE  SITUATIONS 
CINVOLVING  OSCILLATOR  LEVELS  3  OR  HIGHER- 
STOP 
END  IF 
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I F ( CHDENS . LT . STATES ) THEN 

EFERM / =HBAR*CHOENS/RHOarM-SEP/2 . DO 

ELSE 

X 1 =0 . 5D0* ( ( CHDENS/STATES ) - 1 ) 

EFERM1=( 1  DO+Xl )*SEP+SEP/2.D0 
END  IF 

EFEV=EFERMI/CONV 
EFMEV=EFEV* I 000 . DO 
TYPE  8.EFMEV 

FORMAT!/,'  FERMI  ENERGY=- ,F7.3.  niev- ,/) 

DETERMINE  RANGE  OF  ENERGIES  FOR  PL  SPECTRUM 
TYPE  1 1 . EEG 

1  FORMAT!'  EFFECTIVE  ENERGY  GAP  !THEO)=' .F5.3. 'ev' ./) 

TYPE  •.'ENTER  EXPER'MENTAL  EFFECTIVE  GAP  IN  ev ' 

ACCEPT  ’.ECPEXP 

EGPEXP=EGPEXP*CONV 

EM  I N=EGPEXP/SEP- 1 . DO/ 1 0 . DO 

EM AX= I EGPEXP+EBMAX ) /SEP+EFERM 1 /SEP- 1 . DO/2 . D0+ 1 . DO/ 1 0 . DO 
RANGE=EMAX-EMIN 

UPL I M= ! EGPEXP+EFERM I -SEP/2 . DO ) /SEP 

SET  ENERGY  STEPS  IN  PL  SPEC 

STEPS=200 . DO 
ISTEPS=STEPS 
DELTAE=RANGE/STEPS 

DEFINE  PARAMETERS  FOR  INTEGRATION  ROUTINE 

EPSABS=O.DO 
EPSREL=1 .OD-2 
LIMIT=IOO 
LENH=400 

BEGIN  LOOP  TO  PRODUCE  PLSPEC 

NUM=ISTEPS+I 
DO  22  1=1 .NUM 
AI=I-1 

ENERGY! 1 )=EM1N+AI*DELTAE 
X=ENERGY! I ) 

CALL  DQAGS ( F . WM I N , WMAX . EPSABS . EPSREL , RESULT . ABSERR . NEV AL . I ER . 
C  LlMIT.LENIf  .LAST.  IWORK.WORKI 
TYPE  •, ' IER=' . lER. •  ITER=',I 
PLSPEC! I )=RESULT*X*»2 

2  CONTINUE 

DEFINE  SINGLE  PRECISION  VARIABLES  FOR  PLOT  ROUTINE 

DO  75  1=1 ,NUM 
SPLSPEC! I )=PLSPEC! I f 
SENERGY! I )=SEP* ENERGY! I )/CONV 
WRITE!32.» )SENERGY! I ) .SPLSPEC! I ) 

5  CONTINUE 

BROADEN  PL  SPECTRUM  ?? 

TYPE  •,'  DO  YOU  WISH  TO  BROADEN  LEVEL;i  ?  !Y/N)' 

ACCEPT  44, NANS 
44  FORMAT! AI) 
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IF(NANS.EQ. V )THEN 

TYPE  ENTER  BROADENING  WIDTH  IN  mEV  ' 
ACCEPT  •.GAMMA 
GAMMA=GAMMA/l 000 . DO 

CALL  BR0ADEN( SENERGY . SPLSPEC , YB . NUM ) 
CALL  XYPLOT( SENERGY. YB. NUM) 

ELSE 


PLOT  WITHOUT  BROADENING 

CALL  XYPLOT( SENERGY. SPLSPEC. NUM) 

END  IF 

TYPE  •.DO  YOU  WISH  TO  BROADEN  THE  SPECTRUM  ?  (Y/N)’ 
ACCEPT  4 4. NANS 
IF(NANS.EQ. ■ Y' )THEN 

TYPE  •.'  ENTER  BROADENING  WIDTH  IN  mEV  ' 

ACCEPT  •.GAMMA 
GAMMA=GAMMA/1 OOO . DO 

CALL  BROADEN ( SENERGY . SPLSPEC , YB . NUM ) 

CALL  XYPLOT( SENERGY, YB. NUM) 

ELSE 
END  IF 


END  OF  MAIN  PROGRAM 

TYPE  •.'FORTRAN  STOP' 

END 


••••••  BEGIN  SUBROUTINES  •••••• 

REALMS  FUNCTION  F(W) 

IMPLICIT  REALMS  (A-H,0-Z) 

COMMON/CONSTANTS A . A . B . P I , SP I 2 , STWO . SEP , EGPEXP 
COMMON/EXTRA  AJPL I M . GAMMA 

DEFINE  FUNCTIONS 

YY(X.W)=2.D0^(X-0.5D0+EA(W)/SEP)^(A/B)^^2 
AMU(X,W)=DSqRT( 1 .DO+YY(X,W)) 

DETERMINE  IF  Y(X)  >  0  ; IF  YES,  EVALUATE  FIRST  TERM  OF 
INTEGRAND.  IF  NO. SET  INTEGRAND  TO  ZERO  AND  RETURN 


F=0.D0 

Y=YY(X.W) 

tO=YY(UPLIM,W) 

IF((Y.LE.O.DO) .OR. (Y.GE.YO))THEN 

F=0 . DO 

RETURN 

ELSE 

T1=AMU(X.W) 

T2=DEXP( ( (B^Tl )/(STWO^A) )••£) 

GAM=W^B/A 
T3=DEXP(GAM*T1 ) 

X0P=B^T1 /( STWO^A )+W/STWO 
XOM=XOP— STWO^W 

Z1=PI  •B^T2^T3^DERFC(X0P) •  (  ( 1 . D0-(B^T1/A )  •  •2)/2 . DO/Tl  •  •a-GAM/a . DOAi  ) 
Z1=Z1/(2.D0^TI^A) 

Z2=P I  •B^T2^DERFC( XOM )  •  (  ( 1  .  D0-(B^T1/A )  •  •£  )/2 .  DO/Tl  •  •2+GAM/2 . DOAl ) 
Z2=Z2/( 2 . DO^Tl •TS^A ) 

Z3=SPI2^(B^^2)^T2^(T3^DEXP(-X0P^^2)+DEXP(-X0M^^2)A3)/2.D0/(T1^A)^^2 

T4=Z1+Z2+Z3 

F=T4^*2 
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DETERMINE  IF  Y(X-l)  >  0  :  IF  YES  EVALUATE  SECOND  TERM  OF 
INTEGRAND.  IF  NO.  RETURN 

Y=YY(X-1 .DO.W) 

Y1=YY(UPLIM-1  DO.W) 

IF((Y.LE.O.DO) .OR.(Y.GE.Yl))THEN 

RETURN 

ELSE 

T1=AMU(X-1 .DO.W) 

T2=DEXP( ( (B*T1 )/(STWO*A) )**2) 

T3=DEXP(GAM*T1 ) 

XOP=B*Tl/(STWO*A)+W/STWO 

XOM=XOP-STWO*W 

Z1=PI •B*T2*T3*DERFC(X0P)/(A*T1 ) 
Z1=Z1*(B**2*(B*T1/A+W)/(2.D0*A**2)) 
Z2=(PI*B*T2/T3)*DERFC(X0M)/(A*T1) 
Z2=Z2*(B**2*(B*T1/A-W)/(2.D0*A**2)) 

Z3=SPI2*B*T2*(1 .D0-(B*T1/A)**2)/(A*T1**3) 

Z3=Z3* ( T3*DEXP( -XOP*  *2 ) -DEXP( -XOM*  *2 )/T3 ) 

T5=-Z1+Z2-Z3 

F=F+0.5D0*T5**2 

END  IF 

RETURN 

END 

REAL'S  FUNCTION  EA(W) 

IMPLICIT  REAL'8  (A-H.O-Z) 

COMMON/CONSTANTS A . A . B . P I . SP I 2 . STWO . SEP . EGPEXP , EBMAX , SP , ALPHA 
EB(Z)=EBMAX*DEXP(-ALPHA*Z"2) 

WW=W-SP/(2.D0*B) 

EA=EB ( WW ) -EBMAX-EGPEXP+SEP/2 . DO 

RETURN 

END 

SUBROUTINE  BROADEN ( X . Y . YB . NUM ) 

IMPLICIT  REAL'8  (A-H.O-Z) 

REAL*4  X( 1000), Y( 1000), YB( 1000) 

COMMON/CONSTANTS/DX , A . B . P I . SP I 2 . STWO , SEP , EGPEXP . EBMAX , SP . ALPHA 
COMMON/EXTRA  AJPL I M . GAMMA 

LORENTZ(A,B)  =  (GAMMA/PI )/( ( A-B)"2-K;AMMA**2) 

RANGE=X(NUM)-X( 1 ) 

WEIGHT=X(2)-X(1) 

I DEL= 1 1 D I NT( RANGE/GAMMA ) 

ANUM=NUM 

AIDEL=IDEL 

L=I IDINT(ANUM/AIDEL) 

DO  12  1=1 ,NUM 
YB( I )=0.D0 
XC=X( I ) 

J=I+L 

DO  12  M=J-L.J+L 

IF( ( (M-L) .LE.O) .OR. ( (M-L) .GT.NUM) )THEN 

YB( I)=YB( I )+0.D0 

ELSE 

YY=Y(M-L) 

XX=X(M-L) 

YB( I )=YB( I )+YY*WEIGHT*LORENTZ(XX.XC) 

END  IF 

CONTINUE 

RETURN 
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